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1.  Introduction 


riic  rmrcul  iul<.T<'.st  in  liypoi^^ouic  ilo\v.s  invigoratos  tli<’  iiilnixliistipliiiaiy  iv^caich  ijf  chnin- 
ical  kiunticb  aiK.l  acro<iyuaiiii<:.s.  Signilicaiit  progrc^^  lia.s  Ix’cn  nunlc  in  biinulaliny  llin 
non<‘<iuilibiiuni  ^li^^.sociatiun  phonoincnon  lor  higli  tcinpcraluiv  air  as.^'ocialrrl  with  hypor- 
bonir  lliglit  W  lic.'n  lli<'  oncrgy  Icv'-l  ol’  a  ga-soous;  niodiuni  roaclR'.s  an  <'\  rn  lii,uli‘'r  oxcitcd 

htalr.  IIr'  ga:i  aluni.s  ioni/o.  'Iho  moving  cliargorl  parlick':^  ul  tho  ionizorl  ga^  niixUnv.  will 
in  turn,  croato  an  okRlromagnolic  lickl  and  introrliRo  an  artdilional  rliilusion  pro<<‘b^  for 
i.nicrgy  Uanskn"’.  A  r^ocmingly  unrclaUxl  curront  inlorosl  in  rarlar  orobj:;-^;ccUon  rcdiRiion 
ror^oartli  also  r<.'c<'i\<'s  consirkrablo  altonlioik'"'^.  In  tliis  svir’iitilR'  <,lisciplin<'.  a  major  por¬ 
tion  of  the  <'lforl  must  vonccnlralc  on  analyzing  the  scalloring  or  rolk’cling  of  wavo  in  an 
<'l<x  troniagn<.'tic  liold'''^^'‘k  In  this  s<'ns<.‘.  Ur'  two  entirely  rlilferent  areas  of  research  have  a 
common  int<.'rest.  namely  to  ac<[uire  a  better  analytical  capability  for  the  electromagnetic 
liehl. 

Both  th<.’  forc<xl  dilfusion  of  an  ionizcxl  het<.'rogeneous  gasous  jiiixture  and  the  scattering 
and  emission  of  electromagnetic  waves  from  radiating  borlies  can  be  described  b\'  the  tinio 
dependent  Maxwell  expiations.  The  system  of  ])artial  dill’erential  expiations  is  hyperbolic 
and  is  an  initial- value  proljlenT‘.  Since  the  eigenvalues  of  the  system  of  ecpiations  are  ix.'al. 
the  values  along  two  intersecting  characteristics  determine  the  solution  everywheixx  ainl 
any  rliscontinuities  of  a  given  solution  would  be  continued  into  the  interior  domain.  I'or 
all  linite-dilference  time-domain  approximations  of  the  ^.laxwell  equations,  a  fundamental 
dilemma  arises  from  the  necessity  to  impo.se  boundary  conrlitions  on  a  linite  s[)atial  domain 
to  an  initial- value  problem.  In  addition  to  the  possible  degradation  of  the  numerical  solution 
accuracy  of  the  interacting  incident  an<l  scattering  wave  patterns,  the  reflecting  waves  from 
the  artilicial  boumlaries  lead  to  erroneous  accumulations  of  energy  in  the  comiiutational 
rlomain.  leading,  in  turn,  to  unrealistic  modulations  of  the  wave  ani])litude~'^k  Numerous 
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appruxiiiiatioujH  over  the  years  have  Ireeii  <,le\'<'lop<'<l  to  overeoiiie  llii''  rlileimna*'’'* lJi<’ 
most  he<pieutly  adopted  iiuiiierieal  proeedure  is  prohalrly  that  ol  [hKiuist  aiul  Majda*  lu 
their  ap[)roaeh.  tli<.'  uul.)oun<le<l  .surroimdiuu  at  tie-  artilielal  bouudaries  is  simulated  lyv  usiiiu 
alrsoi'hiny  liouarlary  eouditioiis.  In  spile  ot  its  sueeesslul  aj^pheatiojis.  the  api)roxijii<it<' 
l)oLiiidary  eoinlitioii.  ho\\<'\<'r.  still  is  inhen-ntly  liiiiitrsl  as  applied  to  the  total  iK'ld* 
Anotla'r  paeiu.u  it(.'iu  arisiiiy  Iroui  uuuK'rieal  simulations  ol  the  rarlar  si.irnatur*'  analysis 
is  an  uruently  ne<'de<l  inipro\<'nienl  in  nunierieid  <'lhei<‘n<  \''"  ’''''  **’’.  I'or  wane  jjropa.iiittion 
l)heuomenon.  the  a<.le<puite  nunn'rieal  resolution  ol  \va\e  j)aekets  at  a  ^iven  lre<pieney  is 
rlietated  Irv  the  niininial  wave  number  within  ilie  lre«pieney  range.  At  an  extremely  higli 
lre<[Ueney.  e.g..  gigahert.^:  range,  the  reipiirerl  number  ol  rli.seretizerl  data  nodes  to  sohe  a 
l>raetieal  enginwring  prol.ilem  is  enormou.s''-' .  I'he  aggregaterl  eonse<juence  is  that  large 
tunounts  ol'  data  are  to  be  proeessed.  usually  by  a  eonditionally  stable  nuinerieal  algorithm, 
which  limits  our  capability  in  this  area  ol  seieiitilic  endeavor.  .V  jmssible  alternative  nia\'  be 
•  h'riced  from  the  recent  advances  in  llux-spliUing  schemes  lor  solvijig  the  liuler  e<iuations 
in  computational  Iluid  dynamics 

I  he  luudauieutal  idea  ol  llu.x-splittiug  methods  iii  solving  hyi)erbolic  systems  of  eciualions 
is  based  on  the  eigenvalue  analy.sis*'-’'^''^^.  In  linite-clillerence  aj)])roximations.  the  well- 
posedness  re<[uirement  and  th<'  numerical  stability  of  the  solving  scheme  aie  ultimately 
linked  to  eigen vahu’s  ol  the  governing  e<|uatioiis"‘^'^k  I'lie  most  recent  jnogress  in  char¬ 
acteristic  based  and  the  total  variation  dimiiiishing  (T\  D)  schemes^ demonstrate  that 
numerical  stability  and  accuracy  can  be  drastically  im[)rovecl  by  using  an  api)roi)riate  linite 
dill’erencing  lor  sj)lit  ilux  vectors  according  to  the  signs  ol  the  eigenvalues.  In  essence,  the 
system  ol  e<iuations  is  manipulated  to  achieve  the  Kiemann  problem^*^.  hor  the  electro¬ 
magnetic  held  ol  radar  signature  analysis,  the  main  concern  is  the  a])pro])riate  treatment 
of  incident  and  scattering  wave  propagation  from  a  rellecting  body.  Thus,  the  elimination 
ol  nuni)hysical  rellecting  waves  from  artificial  outer  boundaries  is  a  ijaramount  concern 
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('uukl  b<_'  liy  llu;  theory  ul  rluu'aaU.'ii.slic.''' lli<’  luocnl  clluii  w  ill  <i<l<lr<'» 

v'ulh  ksu<'  Ia  u^illy  a  ih-w  K  <i<  ([iur<.'<l  (  'I'  D  t<'<  hui<iii<'  'or  <  li  <  t  i<.aiia,iiii*.'ti<  < dli  iilai  I'Mi:-. 
lilt;  a[)[)rua<;li  lu  uiiiut'rit  al  <;lli<;i<;iK'y  iiu[)ro\fiiifiii  i.>  also  oii  the  uiifouiilftl  t'iut'u- 

wt't.'lur  s^truciiirt'  aii<l  tli<'  iiiUutlurlion  ol  .suilalilt;  iiujilirit  <llyulitllm.'^  lo  i^t.il'f  ilic  Ma.wsfll 
<;t[uali<jii^. 

1  h<;  llux-h|)lluiiiy  liowmt'r.  iiab  an  limilalittu  in  llnil  llif  ttMlln  n  nt 

iiialrict;:;'  ul  tin;  yu\’(;rniny  ('<[ualiuns^  when  wrillen  in  llnx  \<'<.  tt;r  lunii  t  an  !«'  tli.i'^onali/t't.i 
in  tjiily  uiic  <linit;nsiuii  al  a  litiit'*''''"*'^.  Hi<;r«;lur<;.  in  llmltillilIl<'n^ional  analysis.  Uie  ('x<ii  i 
nu-r<;ll<xlion  larlichl  wa\<;  <  un<liliun  i^-  a<;lii<;\'al.)l<;  only  lor  \va\<;  iiioliun  which  ''<;>>(>  a 
duiuinanl  urientaliun.  In  princiiilc.  thit-  lavurahlc  cuntliliun  i:^  atlaiiuihlc  liy  ca.''lin.u  ihc 
Maxwell  e<[ualiunb  in  i;<;neral  curvilinear  coord inateh-.  On  ihe  artilicial  larlield  boirndar\. 
the  courdinale  will  be  adaplcxl  to  ali.yn  with  the  princijial  axi.s  ul  the  wa\e  prupa.ualiun  tt) 
t.'liiiiinate  ur  al  leah>t  lu  uiiniiiiize  unde.sirablc  wave  rt'llecliun.  On  the  interlace  ul  tlit.’leclric 
nic'tlia  ul  dillereiil  i’<;ruiitti\'itieb.  the  cumpk'x  .•scalterin.i;;  ^hajx;  will  l.>e  ])re:;cribe<l  I.a  a  btxK' 
cunlunual  cuurdinate  .surface.  The  .surl'acc  uutwartl  nurinal.  re<p.iired  l.<y  the  !.iuun<lar\' 
rc'laliunshi[)s  of  the  elect ruiuagnclic  lields^h  can  be  easil>’  cuiiijiulcxl. 

In  summary,  tin;  present  iinesligaliun  attempts  to  rlevelup  ellicienl  numerical  pruce<lures  lur 
sul\  ing  the  time-<.l<;pen<lenl  Maxwell  <;<[uatiuns  in  free  space.  Ihe  electric  llux  <.l<;nsity  and 
the  magnetic  llux  density  will  lx;  Iransl'urmed  lyv  a  <liagunalizing  matrix  lu  ac<[uire  invariants 
along  the  respecli\e characl<;ristics.  Ihe  split  an<l  uncuu])le<l  llu.x  \  eclurs  will  then  be  sulvcxl 
by  secund-ur<.l<;r  lemiiural  and  spatially  windward  linite-<lillerence  approximations. 
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11k'  c<rlumn^  of  >  are  dimply  the  eigenvector!?  ol  the  coellici<;nt  matrix  .1. 

IIk'  mullidimein'ional  Max\v<’ll  e<iuationt>  are  i?ei)arate<l  into  one-rlimen^-ional  comiioiK'Hlh 
in  <;ach  <lirection  ol  tin'  ( 'arlei^ian  Irame.  h'or  the  prer>ent  invebli, gallon  ainl  without  lo,'?:-  of 
giaK'ralily.  the  r'Caltered  \va\ei'  are  coii!?i<l<a’e<l  to  he  conlinerl  in  ii?otropic  media  ir^ejiaratixl  In 
j)h\>;ical  interlacia  1.  nd<a'  these  conditions,  tlie  iiermittivPy  and  permeability  are  <issum<'<l 
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loiilliiiU  <'<iual lull.'  !U<hixit<'  lliat  alauiu  I ra |<'<  i<.)ii<.>  wilJi  .''lujx'.s  ul  aii<l  iicUiitu*' 

^^vu<•<l  ul  UuUl.  auvl  willi  Zulu  uiyuiivalviu.  lliu  uhaia<;t<'ii^liu  vai'iabl<.'!^  aiu  imaiuuil.  Mu.'-l 
ini[iui  I aiil ly  lui  llx.'  ,'^vxl<■Ul  ul  <x[iuttiuu;>-  <ill  ll■all^•lul'lll<■<l  ib’ix’iRluiil  \'aiiablu;-  aiu  uuinpluluh' 
I  111'  uu[)l(;<  1  <11 H I  call  l.x;  >ul\  <x  [  iii<.li\  i<lu<i]l\’.  lu  ul  Ix-i  wuiiLs.  t  bu'  iii<ili'ix  .sx  .sU'ui  U'  <iucumpu.''u<  I 
iiilu  u|)[x'i  aiul  luwui  lii<liayuu<il  sti  uui  lll'«'^.  (_)nl\  <i  xin.ylc  su<'<'[)  will  lx;  R'<iuin'(l  lo  llx; 

<  uiu|ilu|u  .''W'lt  .11  ul  <.li;.<.r<;liz<'<.l  u<jiial  iuiia.  l  ui  an  iiiiplicil  ^-uh  iii"  lli<;  ix;ula<ii<i,yuua' 

luvurMuii  pK.'i.ufliiru  Ix.'uuiiiu^  uauuuu.-'aaty  <iii<l  l<;a<l.^  tu  a  <;llR'i<;ut  uuui<;i'i<.al  picx.ciluK'. 
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biiailai  uxpK'^^iuu;;  lui  lliu  Iwu  uiIrt  <;unlpullcul^^  aru  IouikI  in  liie  iilunlical  mauiKT.  1  In; 
uiultiiliuu'iisiuual  p'ubk'iu  svili  bu  isulvub  by  uillxT  a  \viu<l\\ai<,l  uxplkil  lucUiod  or  au  ap- 
piuxiiualu  iauluml  implicit  :sclicmc'  '"'x  A  spccilic  isclcctiun  ul  numerical  algorilluiLs  is  one 
ol  the  major  reM'arch  goal:;,  criterion  ol  the  choice  i:;  that  the  algorithm  yield  low  <,li:;bi- 
[)ati\e  and  <.li8[x'r,si\e  error.':'.  In  ihi:^  connection,  mmierical  .'?cheme.s  which  can  .■^ati.'^ly  the 
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i'oii<,lit ion  ill  tlio  <l»'!|<.'ii<'i'al«.'<l  oii<'-<liiii<‘ii^iioiial  .■'<-ttiiig  is  suiiulil"  ili<'  case  vl 

iiuiik'iiK'iilaliuii  ol  IjouiKUiry  <.uii<lilioiis  uii  rillR-r  a  s<;atl<'i'iii,y  body  oi  llio  ouloi  <'oiiijHit<i- 
lioiial  (^loiiiaiii  is  also  a  iiuijor  (.oiisidoiatioii.  l-'or  a  siiiijilo  waw  oqiiatioii.  tlio  [iicloiriK  <• 
s<'<Miis  to  lia\<'  l.xK'ii  l.>ias<'<l  toward  t!i<’ explicit  wiiulward  iiietliod  Iry  W aniiiii,!;  and  13<’<iiu' ’ 
llow<'vei.  <1  wiiidwairl  implicit  iiK'lliod  also  possoses  this  niiiquo  pKjpt'ity  and  is  iiiclu<l<'d 
ill  the  present  in\ e^^t iyat ioir'^. 

A  critical  ste[)  li<'r<'  is  that  the  liiiit<'  dilh-reiiciiiy  is  iinpleineiiterl  to  honor  the  physical 
ori<‘ntation  of  the  wa\<'  propayat ioir  ”’  *'^.  hor  example,  lorwarxl  aiirl  backward  (.lilh-reiicing 
of  the  spatial  \arial.>les  are  lornied  lor  the  iieyalive  and  the  positi\’e  eiyenvahie  respectively. 
Ill  ,ii,<'ii<'ral.  iiilornuUioii  is  [)ropagal<'<l  I'rom  the  interior/exterior  ol'  the  coiniwitalioiud  <,lo- 
iiiain  to  th<'  boniidary  akiiig  <haracteri.stics‘  *x  Ihe  characteristic  e<iuatioiis.  l:i<i  (b).  are 
r<'all\'  th<'  coni|)c.til)ility  r<'latioiiships  ami  are  <.'<iually  ajiplictibh’  on  the  boumlaries  of  the 
coinpntat  ioiial  (loinaiii.  lo  use  characteristic  data  to  establish  the  exit  wave  boum.larv  con¬ 
ditions.  knowledge  ol  the  orientation  ol  the  wave  I'ront  with  resjiect  to  tlx.'  boumhiry  surface 
becoiiK.'s  essential.  .V  i.nii<[U<'  leaturi'  of  tlx' electroiiiagnetic  Held  provides  this  re<|uired  in- 
forniation  nicely.  In  free  space,  tlx-  direction  of  wave  jiropagation  is  always  known  since 
tlx'  eleclroinaguc'ti(.'  waves  move'  in  tlie  direction  perjiendicuiar  to  both  the  electric  and 
magnetic  lields.  L  •  //'  ‘b  .Near  the  outer  control  surface,  the  characteristic  e<iuations  will 
lx.'  project<.'d  to  the  j)riiicipal  axis  of  the  wave  inolions  based  on  the  known  local  solution. 

covu'dinate  translorniation  Iroiii  the  Cartesian  Iraiiie  to  a  general  curvilinear  system  can 
then  be  a<lopted  to  ensure  coordinate  alignment  with  wave  motion,  .\nother  advantage  of 
the  coordinate  transforiuation  is  that  it  permits  a  high  numerical  resolution  by  local  mesh 
n.'lineiiK'iit  if  nec<'ssary'b  .Nkiatiwhile.  the  compk.'x  scatt<.'ring  body  is  easily  represented  by 
l)ody  coiiformal  coordinates. 

general  curv  ilinear  coordinate  transformation  is  introduced  by  delining  a  one-to-one  rela¬ 
tionship  between  two  sets  of  independent  variables.  In  the  present  analysis,  the  coordinate 
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.s}  liiiiiictl  lu  tlic  spaluilly  \aiiablo. 


^  -  a.f. ,//.:) 
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A  unc-lu-uiK'  <uor<jiiiat<'  Uatisluiuiatioii  ciimikxI  ihiuuyh  Ur'  iiunvaiiij^lini"  Jacobian  ol 
<our<liuate  ll■all^^lullualiou.  Uy  iIr'  chain  rule  ol  rlillcrcnlialion.  iIr'  Maxwell  cxjualionh. 
ac<juire  the  lollowiny,  loriu.  similar  to  that  in  the  (.’arte.sian  .system. 
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Uk'  tluee-<,limensional  proljlem  i.s  separated  into  three  one-rlimensional  problems  along 
<'ach  traJisloriiK'd  cooixlinate  lin<-.  I'or  the  present  purpose,  only  one  comjwnent  ot  the 
Maxwell  <'(iuations  is  prt.'sented.  ilm  other  components  are  nearly  identical,  dillering  only 
in  the  c\clic  rotation  of  th<’  metrics  of  coordinate  transformation.  The  eigenvalues  ol  the 
geiu.'ial  cuiR'ilinear  .systi.mi  are  obtained  after  some  algebraic  manipulation.  As  expected, 
the  functioiial  forms  are  similar  for  all  components  in  the  transformed  space  and  the  details 
will  not  be  pres<'nted  h<.'re. 
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I'ullowiug  a  similar  plot  r.luiv  in  do v< -loping  llio  oliaracloii^tio  otpialions  loi  llio 


( 'atlt'^iaii  a  lypn.al  .-^ol  ol  iiuaiianl^'  along  oiio  <  ooi<linalo  i^. 
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llio  nuniorioal  solution  of  llio  Ma.wvoll  o<ination  it>  obtainod  in  oliaraotorisilic  \ariablos  by 
lliroo  dilloronoing  opoiatons  along  lian.sl'onnod  coordinato  linos,  r'orwanl  or  backward  linilo- 
tliiroronco  a])pro.\iniations  aro  lorinod  according  to  llio  sign  ol  llio  associalod  oigoiualuo. 

M  llio  inosonl  stage-  ol  algoritliiii  eloxolopnionl.  only  two-diinonsional  siinnlalions  aro  al- 
U.-niptod  lo  so<‘k  conlirnialion  ol  llio  basic  concopt.  Tosl  casos  consist  ol  incielont  trans\orso 
oloclric  aiK.l  niagnolic  wa\os  proiiagaling  llirougli  a  linito  conipulational  domain  ami  oloc- 
ironiagnotic  pulses  gonoralod  by  an  oloctrical  curront  .source-  j>e-rj)e-nelicular  to  llio  plane-  ol 
wa\<-  motions.  .Mlliougli  llie-.so  simulate-el  phe-nome-na  are-  e-lomonlary.  tlioy  pose-  soeoroly  ele-- 
maneling  numorical  accuracy  re-fpiiromonls  lor  ele-scribing  eliscontinous  wa\e'  boliax  ior  anel 
yot  possess  all  ossontial  le-aturos  ol  a  scalle-ring  li<‘ld.  Tlie-re-loro,  the-  physical  boundary 


foinlitioiiis  ul  a  scallcriiig  lju<ly  arc  iiol  altciiiptc*.!  in  the  jncbciit  analybib.  1.  iidcr  tliit) 
rraiiicwork.  luiiucrical  rcjsuUb  rcllcct  ihc  i^csi  pussiilrlc  pciioriuaucc  ol'  the  new  nuiiK'iieal 
prueerlnre^  in  ^iml.llalin,^  tlie  ^^eallerin^  eh'eliuniagnetie  liekl. 
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3.  Nunieiical  Methods 


SiiKO  Maxwell  <.'<|aaliuii,s  aie  <liagoualiz<'<l  iu  cadi  .sjialial  ‘lirccliuii  au<l  cuuij)cl<'lcly 
uiKuupkxl  ill  tcMii.'s  ()1  ci.ycin ccioijs  aiwi  a.s.sodat<xl  cigcii\'aluc,s  (positiw.  iicgali\c  and  null), 
tlic  windward  linilc  (.lillcrcncc  ap[)roxinialion  is  uscxl  to  del  ter  (.leseribe  llie  physics,  llie 
solution  proeerlure  lor  the  ei.geiueetor  ikxl  to  the  null  eigeinalue  (E<[s  (1.  an<l  1(1.17)  is 
trivial,  i.e.  the  unirliixxtion  eigenvector  is  an  invariant  with  respect  to  tiiiK'.  lor  a  two 
dimensional  problem,  only  two  split  llu.x  expiations  with  nonzero  eigenvalue  are'  nee-eleel  to 
l)e  soh'i'el  in  each  coorelinate  elirex  tion.  'Ihree  s<xon<l-order  tenijioral  anel  spatial,  explicit 
anel  implicit  numerical  sche'iiies  are  implementexl  lor  the  juesenl  iinesligalion.  .Ml  the 
algorithms  consieh.'red  are  either  basexl  on  or  relateel  to  Warming  anel  Beams  works  lor 
solving  the  Euler  expiations"’"'^. 

1  he  two-step  uiiwinel  explicit  methoel  by  Warming  anel  Be.'anr’'  has  bc'cn  aelajvtexl  lor  solving 
the'  Maxwell  expiations’.  This  c'xidicit  scheme'  has  a  conditional  stability  proiie-rty  that 
re.'stricts  the  allowabh.'  time  slc'i)  by  a  (.'EL  value  ol  2  lor  a  one-elime'nsional  problem.  The 
theoretical  limit  is  re'elue'<.'el  to  a  value  ol  2/ \/d  in  lhr<.'e'-elinie'nsional  simulations.  In  jiractical 
applications,  the  allowable  time-sle'p  .size  which  is  elictateel  by  the  minimum  griel  s])acing. 
can  be  elrastically  nxluceel  even  more  by  the  thin  film  i)artitions  ol  elielc'ctric  mc'dia. 
dnc'-elime'nsional  solutions  ol  this  particular  scheme  sat  is  ly  the  shilt  conditioiE''^'^'^  at  a  (.'EL 
nunilx'r  expial  to  1  anel  2.  Solutions  that  poss<,'.ss  this  uniepie  projie'rty  contain  no  ejuasi- 
physical  error  and  are  therlore  highly  elc'sirable  lor  accurate  elc'scrijition  ol  rellccte'el  wave's 
in  the  lar  held  ol  a  scattering  problem.  In  the  prc'sent  lormulation.  the  two-ste'])  methoel 
can  lie  simpliiieel  to  a  single-step  algorithm  by  substituting  the  ])re'elicte'el  result  into  the 
corrector.  Die  resultant  finite’-elille.'re'nce  appro.xiniations  lor  the  positive  anel  the  ne'gative 
eigenvalue  expiations  aexpiire  the  lollowiiig  lorms: 
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I'ur  tlic  uuilti<.liiii<.uisioual  probltnu.  llic  cyclic  Iraclioiuil-slc])  se<piencc  and  the  total  sum 
{jicccdure  weie  iuclu<,le<,l  in  the  ])reseut  analysis.  Numerical  solutions  generated  from  both 
procedures  r<.'taiu  second-order  accuracy  in  space  and  lime.  There  is  essentially  no  disi)<u  - 
it\'  between  solutions  by  these  two  e.\])licit  methods,  although  the  fractional-step  method 
retiuires  additional  data  at  tlie  intermediate  teni])orai  step.  In  addition,  the  numerical  elli- 
ciency  advantage  of  the  liaclional-stc})  method  over  the  total  sum  scheme  no  longer  e.xists 
for  sol\  ing  wave  eciuations.  Therefore,  no  further  elfort  was  devoted  to  its  evolution  at  the 
present  time.  For  the  e.X])licit  methods,  the  numerical  boundai'y  conditions  at  outer  edges 
of  the  computational  dourain  are  easily  implemented  by  specifying  data  of  the  exact  in¬ 
coming  wa\  e  compcueut  and  by  enforcing  the  null  value  of  the  eigenvector  in  the  directioii 
of  the  exiting  wave. 

I'dllowing  the  development  of  Waiiuing  and  Beam’s  trapezoidal,  second-order  implicit 
approximale-factored  algorithms  for  solving  the  two-dimensional  hyperbolic  systeni'^^,  the 
linite-rlifference  api)roximation  can  be  given  as  the  following; 
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OiR-  iiulo  tlial  during  dillninit  numerical  ^.wceph.  liu'  rlilleivm  ing  ojx.uaiors  du  nol  a<  I  on 
tin.'  same  <.'ig<’ii\<xtor.  1  horelorc.  alter  <.‘a<;li  nuiiKuieal  rswr-ep.  lli<’  eigemectur.''  a.'^^o<  ial<’<| 
with  the  next  spatial  urienlalion  need  to  lx-  updaterl  l^y  the  niot^l  ree<'nt  rak ulation.  lids  is 
the  direct  conse<iuence  ol  our  inalrility  to  simultaneously  diagonalize  all  coeliicient  matrico 
oi  I  lie  Maxwell  equations.  Since  relationships  linking  tlx'  eigenceclors  ainl  the  electromag- 
nelical  held  variables  are  explicit  an<l  unambiguous,  this  inconvenience  does  not  impose  an 
inhil.iiting  constraint  lor  lurther  de\eloj)ment  into  three  rlimensions.  On  the  other  hainl. 
the  anticipated  I’uture  generalization  to  inclufle  arbitrary  cuiuilinear  coorrlinates  and  di¬ 
electric  merlia  will  incur  a<l<.litional  errors  ol  linearization In  the  pre.sent  in\estigatio]i. 
the  system  ol  e<[uations  is  linear  Ijv  virtue  ol  the  lacts  that  all  phenomena  occur  within  an 
isotropic  iiKxliumand  all  metrics  ol  coordinate  translormation  are  assigned  constant  \alues. 
1  hen.'lore.  th<.'  similar  matrices  ol  diagonalizatiou  are  iiidc.gx.'iident  ol  space  and  time  and 
ran  b<.'  luought  into  the  dillerenlial  operators  to  complete  tlx.'  lormulation.  Otherwi.se.  the 
linal  characteristic  r.'rpiations  can  be  achieved  only  alter  a  local  linearization  j)rocess.  'I’his 
is.  howe\<.'r.  an  inlx.’rent  h-ature  ol  using  implicit  methods  to  sohe  a  nonlinear  ecpiation 
system. 

In  the  original  work  ol  Warming  and  Beam  on  the  construction  ol  imj)licit  lactored  scln.'mes 
lor  hyperbolic  systems''’,  a  trapezoidal  lormula  in  tinx.'  aix.l  internally  inconsistent  wind- 
warrl  scheme  was  establislx.'rl.  Internal  inconsistency  iini)lies  that  the  explicit  and  implicit 
rlillerencing  lor  temporal  evaluation  are  unbalanced.  The  <.'X])licit  (Kllb)  dilierencing  is 
gi\en  l.)y  a  three-point  one-si<led  approximation,  whereas  the  implicit  (LIIS)  dilierencing  is 
re])r<.'sente(l  by  a  two-point  one-sirh.'d  a])pro.xiinatioii.  Although  the  latt<.'r  dilference  aj)])ro.x- 
imation  is  lirst-orrler.  it  operat<’s  on  the  increment  oi  characteristics,  and  as  a  con.se<px.nce 
the  resulting  scheme  is  spatially  second-order  lor  aJiy  li.xed  ratio  ol  temporal  aiul  spatial 
increments-'''.  This  imi)licit  scheme  is  conditionally  stable  lor  the  trapezoidal  lormula  and 


13 


it  :>ali^lics  iIk-  isliill  cuntliliuu  at  (.'I  L  <.-<[iial  tu  1.  IIk-  [uociit  tiuR'-ceuteicd  and  con.sisiicut 
tlu<.'<'-puiut  wimlwaixl  .scla'ino  i.s  (.kxclojxxl  irnni  tlii.s  1  rat n<' work,  dhe  ])oint  ol  <.k‘j)arturo 
li*'^  in  the  tact  that  the  guserniny  hvjx'ibolic  ^)y^ileln  of  exjuations  is;  uncuuphxl  by  diag- 
oniilization.  Ihc  pontatliagonal  coellicicnl  niatiix  struct  tiro  of  a  bidirectional  windward 
rlilierencing  ap])roxiination  is  eliminated,  lu  tlie  present  sj)lit  llu.\  implicit  formulation, 
each  characteristic  cxpialion  cam  be  cast  into  a  simph,'  bidiagonal  matri.K  structure  and  the 
imx.'rsion  process  is  straight lorward.  liie  sp<xilic  expiations  sohed  arc; 

X  sweep:  I'or  A  >  0: 

<(..;■ -2)) +  (.ia-;(ny)- 

+  1)  +  »V(/-2+2))]'‘  I 

1  dr  A  <  0: 


Idr  A  —  0; 

tdi'-A)  = 


(‘^)  I'C*'  +  l-j)  +  <(i  + 

H^)  [«■;.(' +  2.a  +  <(.  +  2,aj- 

j-i)+ 

-2))  +  (■iwp'-ll- 
+  J)  +  +  2))]  ' 

i+d«)) 


(3.6) 


-  1))  +  +  1)  -  ■w~{uj  -  i)] 


(3.7: 


11 


Y  sweep:  I'oi  A  >  U: 
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111  general,  tlie  algorithiii  is  uiicoiRlitioiially  .slaljle  wlieii  applied  lo  a  simple  wave  e<iual.iou. 
a  feature  which  may  be  \ery  \aluable  in  improxing  numerical  elliciency  for  computational 
techni<iues  in  electromagnetics,  'i’lie  elementary  stablity  analysis  of  this  algorithm  is  ])re- 
sentecl  by  the  modulus  of  amplilication  factor  and  the  relatixe  j)hase  error  in  rigures  1 
and  2.  respectixely.  The  xalues  of  moxlulus  are  depictcxl  for  the  range  of  ('FL  numbers 
from  0.7b  to  b.O.  It  is  clear  that  the  numerical  method  is  dissijiatixe  an<l  the  numerical 
error  is  diminished  monotonically  as  (’FL  is  increa.sed.  On  other  liaiul.  the  relatixe  phase 
angle  e.xhibils  a  predominant  leading  error  for  (.'FL  xalues  less  than  unity,  then  sxxitches  to 
a  iiersistent  lagging  error  for  (.'FL  values  greater  than  2. 
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I'igvue  i;  Modulus  of  Amplification  Factor  0.7o  <  CFL  <  5.0 


In  llin  iang<'  uf  ( 'i  L  nuni!joi.s  truni  1  to  2.  iho  rclati\'o  pliat^c  angle  changes  from  a  slight 
leading  error  to  a  lagging  error  aecurding  to  the  increasing  wave  luiinber.  lire  eross-cner 
ol  rc-lativc'  phase  angl*.'  from  leading  to  lagging  c'lror  takes  ])laee  in  the  mid-range  ol  wa\'e 
numlxn  and  constrains  the  discrepancy  to  a  relative  low  magnitude'  in  a  wide  spectrum. 
•Most  importantly,  the  re'lative  i)ha.se  error  is  signilieantly  less  tlum  most  time-  and  space- 
centered  iniiilieit  schemes  in  the  high  wave  numljc'r  clomain'^"*^. 

for  the  appro.ximate  laetored  implicit  .scheme,  the  implementation  ol  numerical  boundary 
conditions  at  tlie  intc'rmediate  time  step  is  more  eomjdiealed  than  for  the  ..\plieit  methoct. 
In  til*'  ca.s<.'  when  an  e.\act  datum  at  the  artiliehil  boundary  is  known,  the  intermediate 
value  must  Ice  <'.\trac  t<'d  from  the  ne.\t  consecutive  numerical  swec]).  That  i:..  the  boundary 
conditions  of  e<[uations  CJ'J).  C-’d).  and  (21)  must  be  derived  from  e<iuations  (2b).  (2b).  and 
(27)  respectively,  wliereas  tor  the  e.-cplicit  calculations,  this  issue  is  completely  eliminated 
i.)y  t!ie  total  sum  procedure,  kor  the  ease  where  a  nonielleeting  wa\e  condition  is  recpiired. 
tin.'  null  value'  of  the  c.'ige'ii vec  tor  opposite  the  e.’citing  wave  is  imjeosed.  Them  any  error  that 
may  b<'  incurred  is  contrilnited  Icy  the  misalignment  of  the  orientation  of  wave.'  propagation 
and  the'  coordinate  direction.  In  principle,  this  error  could  be  remedied  by  the  coordinate 
transformation. 
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4.  Discussion  of  Nvimerical  Results 


All  nuiiR'ikal  results  were  •i.oucralel  on  a  la  •  la  vsjually-spacut]  Cartusiaii  s\  sU'in  in  wlildi 
iiu'l  iias  ol  iho  <;oor<iiual »'  Iraiislornialioii  assuiucrl  a  \<tlu<’  ol  <’ilh<‘r  1  or  <1.  *  has  ili<'  liii- 

cari/atioii  error  is  coiiiiilelely  liniinaterl  Ireun  the  pnx'iii  aiial> sis.  l  urtiK'r  «i«'\<'lu|)un'iii  <>1 
I  Ik'  present  iiiei liorlology  lor  praet  ieal  <ippliealions  to  iu<  lu<le  <|i<  l<’eiri<;  iiHslut  <iii<l  arl>itiar_\ 
scat t (.'rill, •;  hurlnrs  will  lie  [mrsiusl  in  upe.juiiny  <-llorls.  An  lUlS  llJ/i'Jlt  (11  \  workslaiioii 
was  useij  exelus i vely  lo  generate  all  muiieri<  al  results,  ll'e  relative  uuuierieal  ellieiiuiey 
of  all  stiK.lierl  nunierie  'l  proee<lur<'s  is  proUibly  lier'  evaluated  by  the  eount  ol  aritluiK'tie 
operations  to  proeess  the  two-dinn  nsional  liehl  <lata. 

|t  is  olivious  that  the  upwind  e.'iplieit  m<'tho<l  re<iuir<'s  tin*  hwvest  arithiiielieal  opeiiitions 
to  eoiuph'ti' a  lie*'l  point  ealetilaiion.  Howev<‘r.  th<'  trapezoidal,  eonsistent  iinpli<  it  iiietho'l 
has  no  iiu|)o.sed  reslrietioii  Iroiii  the  staliilit.v  eoiulition  and  can  juoeess  data  at  far  .yn.'ati’r 
time  steps.  I  his  lavoralile  nuiin'rieal  leature  will  he  <’ver  inure  apiiarent  wlien  the  <liel<'etrie 
nnslia  ar<'  sei>arate<'  lyv  a  thin  eoatin.!;  an<l  tin'  time  step  tilluweil  is  eont  rolled  hv  the  smallest 
spatial  vlimension.  Supported  hv  a  suitaiile  mesh  system  for  requirevl  numerieal  rrsolution. 
the  trapezoidal  eonsist<'nt  Implicit  m<’(ho<l  eoukl  have  a  mueh  higher  numerie<d  ellieieney 
than  most  proeevlures  eurri'utly  in  us<'''  ' 

fable  1:  .\rithm«'tieal  Opi.natioiial  Count  ol  Aunierieal  Sehemes 


fjasie  algorithm 

.\rithmeiieal 

Operational  eount s 

Single-step  I  ]nvind.  E.xjilieit 

■yj 

irapezoirlal.  Ineonsistenl  liui>lieit 

76 

Irapezoidal.  Consistent  liuplicil 

io;3 

If) 


.\ii!U(;iical  ol  ^I'uujf.s.  1' ii>l  a  .Muicb  ol  uii**  <.liiiL<.'u;!>ioual  \va\'<.‘  uiotiaiu 

Citk'ulatioiis  t(,u'  all  tlu'ac  miUK.'i'K'al  lio<.k  iU  a  i'aii,‘;<'  ul  (,'1'L  x'aluos  Iroiii  U.Ta  to  ll.U 
i>  .\ii!ii<'iic<il  o'.s’ilt.s  ol  explicit  and  implicit  mc'tlio<is  an.'  imc.stigatcd  for  tlicir 

ali'lity  to  a<.(  oiiioiUitc  the  pio|><iyat  ion  ol  <»  pi<'<.<.'\vi.'?<‘ contiiiou.s  dastuiljancc  ami  to<|aautily 
<["<i>i-[)liysi(.al  cr[oi>.  In  th<'  .''<'<  <.in<l  yroip)  ol  n-Milt^.  tl■au^^\■c^s<■  wa'.c^;  tia\cling  obli<p!cly 
across  ortho.aoiutl  two  <.lilU('ll^ional  cour<linal<'>  an'  pn.',s<'nt<.'tl,  ili<.'  iiuiiiciical  !^cil<.'lla.■s  <in.' 
<'\ainin<.'<l  lor  ^enr'i!  i\  it  \  ol  ccauputc'<l  wa\«'  struct  nn-  n'l<iti\'<.'  to  tin.'  impo.scrl  bt-nmUiiy  coii- 
<lilion8.  i'  inally.  the  new  iiiiuieri*  <tl  m<.'tiio<.l.':.  an'  .ippli<.'<l  to  ^inlnlate  a  locally  nonaiialylical 
oniwan.l  pitjijayaliny  elect  roui<i,yuetic  liekl  .iyeueiated  by  au  electric  current  j)c'rpc'uclicular 
to  the  plane  ol  tran;'\<'i>e  wace^.  In  t hi.''  ,!;ronp  ol  re^>ult^.  the  conil^inerl  ell’ccts  ol  a  sin¬ 
gular  pert  It  illation  ami  its  outward  propa.tiatioii  throuu,h  a  linite  compuational  domain  are 
di'liiK  'alerl. 

In  1'  itrun.’  -J.  tin'  <'xact  <.'lect  rical  In'ld  density  ol  a  traveliuj;,  wav<t  is  comparcel  svith  nunicrical 
results  Iroin  the  traiiezoirlal  consistent  iniiilicit  (T(.'l)  ami  the  single-slc’j)  ujiwiml  exjilicit 
(.SI  11 )  schemes  at  the  (.'I'L  \aluc  oil’.  Ihis  is  also  the  maximum  allowable  time  step  that 
cati  be  specilied  lor  the  SI  E  methorl.  1  he  .solutions  are  j)re.sentcn!  at  itistances  when  a  right- 
runnmg  wave  Iront  reaches  the  mid-iioiut  ol  the  computational  domain  (Fig.  da)  and  exits 
the  numerical  boundary  respectively  (Fig  db).  I'iie  shift  i)ro])erty  of  the  solution  by  the 
SI  L  schenn.'  inrlicates  a  peiiect  translation  of  the  initial  value  in  .sjiace.  The  rliscontinuous 
incoming  wa\e  Iront  is  captun.’d  precisely  by  the  SEE  method  but  not  by  the  'K.'l  scln'me. 
When  the  imiiulsi-  wave  front  moves  through  the  computational  domain,  the  SEE  method 
rluplicates  the  exact  solution  at  each  am!  ewery  rliscretized  point  (Fig  da).  This  highly 
rlesirable  j)ro])erty  of  a  numerical  solution  is  not  preserverl  for  general  two-dimensional 
computations.  In  contrast,  the  solution  of  T(.'l  .scheme  yielded  au  S.l-percent  maximum 
rliscreiiancy  from  the  e.xact  solution. 

In  spite  of  the  fact  that  the  TEl  method  has  no  theoretical  stability  limit  on  the  lime 
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llic  iiuiiK'iical  accLinKy  cuulrulWl  by  tlio  lim<-  step  adopted,  lii  tlio  present  inves- 
ligalioti,  the  <lesire<l  nuiiierieal  aeeur^iey  ean  Ire  iiiipru\ed  by  ledueing  the  (.'I’L  number, 
ur  e<iui\alently  by  adrling  gri<l  norles  to  r<-.sol\e  a  given  wawleiigtli  while  keejring  tlie  Cl'h 
\alue  unalt<'re<.l.  In  I'igure  1.  the  enrielnxl  spatial  resolution  is  achieved  by  changing  the 
periorl  ol  the  wave  Irom  Irr  to  'J.t.  instead  of  rloubling  the  number  ol  nodes  in  space  from 
1')  to  t>(l.  As  a  consecpience.  e<ich  wave  will  l»e  r<-solv<xl  by  hr  points  in  comparison  to  the 
pri'vious  d'J.  <uid  the  ma.\inium  numerical  error  at  the  same  ( 'FL  value  ol  1  diminishes  to  a 
magnitiKle  ol  U.IS  percent.  In  essence,  if  the  Ih’l  .scln'me  has  a  suitable  sup])orting  mesh, 
a  great<'r  Ilexilrility  exists  in  temporal  evolution  of  accurate  numerical  results.  By  virtue 
of  the  mor<'  favorable  stalrility  property  than  explicit  schemes,  the  allowable  lime  sle])  ol 
the  imt)licit  schem<'  ch'peuds  only  weakly  on  the  mesh  spacing.  'Ihis  particular  properly 
permits  a  wid<.'r  time-ste[)  selection  and  jrerhaps  will  lead  to  a  more  ellicient  procedure  lor 
problem  solving. 

In  Figure  O.  the  exact  solution  of  the  electrical  held  intensity  and  solutions  irom  the  SUE 
and  K.'l  sclK'nies  are  presented  together  for  a  (,'FL  of  1.-3.  At  this  condition,  the  ])erl'ect  shift 
doc's  not  prevail  for  the  solution  via  the  .SI,  F  method  at  the  discontinous  wave  front  (Fig  3a). 
Numerical  solutions  from  th<'  two  methods  indicate  distinctive  behavior  at  the  singularity. 
Ihe  solution  via  the  T(.'l  method  approaches  the  undisturbed  held  monotonically  while 
the  solution  v  ia  the  Sf  F  method  yields  an  oscillatory  j)altern  near  the  linite  jump.  As  the 
wave  front  passes  through  the  numerical  domain,  the  solution  by  the  SFF  method  indicates 
a  suj)erior  numerical  accuracy  over  the  .solution  by  the  T(.’l  Scheme.  The  deviation  from 
the  exact  result  is  less  than  ‘i-percent  (Fig  3b). 

The  comparison  of  solutions  of  SFF  and  T(.'I  scheiin's  and  the  exact  solution  for  a  simple 
wave  with  a  period  of  2-  and  CFL  of  2/  \/3  is  jvresented  in  Figure  6.  'Ihe  selected  (.'FL  value 
coincides  with  the  theoretical  limit  of  the  SUE  method  in  three  dimensions.  Again,  numeri¬ 
cal  errors  at  the  singular  moving  wave  front  are  observed,  but  the  numerical  behavior  of  the 


Iwu  m<.'lho(,l!:i  i.s  nncisod.  Now  tlic  SI  E  ^solulioll  appioaclios  tlic  uiKlissliirbotl  vkrliical  lickl 
uioiiotoiiically  ami  iIk-  iCI  .'^okilioii  oxliibiUs  oscillalioiis  across  the  jump  with  a  maximum 
noniializol  discrepaiKy  ol  b.b  percent  (Ei“  tia).  In  the  entire  domain  downstream  ot  the 
wave  Iront.  the  numerical  errors  ol  both  s<  hemes  are  conlined  to  a  magnitude  arouiul  a 
one-luuKlrerllh  ol  oia;  [xiceiil  ( 1' ig  (ib).  I  roiu  these  <)bser\;\ii<.u!'-,  i  lie  optimum  application 
range  ol  the  ICl  sclienu'  lor  high  numerical  resolution  electromagnetic  simulation  may 
resi<le  within  the  ('kh  \alnes  al,)o\e  unity. 

In  Eigure  7,  lor  the  lirsl  time,  all  numerical  .solutions  develoj)ed  in  the  present  analysis 
are  compan-d  with  the  exact  .solution  at  a  CEL  ol  unity.  E'or  the  trapezoidal  inconsistent 
implicit  (  ill)  .scheme,  this  CEL  \'alue  represents  the  maximum  time  slej)  allowable  from 
the  stability  analysis  ol  a  linear  initial  value  system.  Both  solutions  ol  the  SEE  and  Til 
method  poss<'ss  the  shill  i)roi)erty  and  thus  j)roduce  no  numerical  error  at  the  mov  ing  wave' 
Iront  (i'ig  7a).  it  may  be  ol  interest  to  note  that  the  Til  schejiie  is  oid\'  oiu.’  ol  a  lew 
implicit  algorithms  that  can  produce  the  sim])le  wave  solution  with  the  shilt  properly.  The 
solution  ol  'E(,'l  scheme  alone  shows  numerical  overshoot  j)receding  and  dissipative  error 
trailing  the  moving  wave  Iront.  .\s  anticipated,  solutions  of  the  SEE  and  the  Til  schemes 
reached  perlecl  accord  with  the  exact  solution  in  the  complete  solution  domain  (Eig  7b). 
IIk,'  solution  ol  'ECl  scheme  uiKlerpredicts  the  am])litude  ol  wav'e  by  a  maximum  of  1.7 
percent  and  contains  an  identihable  leading  phase  error  which  has  been  verified  by  our 
('arlh.'r  discussion  ol  E  igure  2. 

'Ehe  hist  comparison  ol  all  numerical  solutions  is  depicted  in  E’igure  8,  where  the  CEL 
number  is  asigned  a  value  ol  0.7-3.  At  the  discontinuous  wave  front,  all  numerical  simulations 
exhibit  oscillatory  error  across  the  jump.  The  numerical  errors  normalized  by  the  exact 
wave  amplitude  span  an  extreme  magnitude  ol  i0.*2,  10.0,  and  13. E  percent  respectively  lor 
the  SEE.  Til  and  TCI  schemes  (E’ig  8a).  In  the  postwave  front  region,  the  solutions  by  the 
SEE  and  Til  schemes  are  nearly  identical  (EEg  8a).  The  largest  error  is  again  exhibited  by 


lUc  K.'l  scIk.'uk-'  auU  aiiiouiilii  lu  ihau  2  poixanil.  llowcvt'i.  one  ssliouW  iioU.'  that  lia' 
K.'l  .srla.’iiio  if?  not  i(.‘<,x)niiiR-n<i<.'<.l  iu  Uii.s  of  ( 'l-'L  values  in  practic  al  a]>pii<  aliun.  Hie 

newly  <,le\elope(l  trapezoi<lal  eunsislent  iinplieit  proeerlure  may  be;  eritieal  lor  iniprosemenl 
ol  nunierieal  ellieieney  of  eleelromagnetics  simulation.  Overall,  the  ICl  selieme  is  easily 
the  Ijc.'st  ehoiee  lor  lurther  rlexelopment. 

In  the  secoinl  group  ol  nunierieal  re.sults,  the  three  nuiiK.'rieal  proeerlures  are  eomjiarcxl 
lor  a  j)ro])agating  trans\erse  \va\<.'  at  an  angle  ol  iiieirlenee  to  the  eoorelinates.  In  all  eases 
sturlieel.  the  e.xaet  solution  ol  the  ineirlent  wave  was  u.seel  to  o\erspeeily  nunierieal  bounelary 
eomlitions  ol  the  initial  \alue  system.  l.’iKler  this  eoiKlition.  tin.'  we'll-iioseelness  issue  ma> 
arise'  lor  the  (.lilferenee  system""*,  'ihe  prescrilx'el  exaet  wave  solution  at  the  truneaterl 
Ijoundaries  may  still  appear  as  a  perturbation  to  the  elillerenee  system.  The  jiossilrility  ol 
sustaining  a  eontinuous  nunierieal  solution  by  a  wklcT  class  ol  ])erturbations  to  the  initial 
data  was  assessed  here.  All  .solutions  were  evohed  lor  more  than  a  thousand  time  stejis  at 
a  l.TL  ol  unity  which  eorresponds  to  more  than  22  cycles  ol  a  right-running  wa\e.  The 
implicit  schemes  ar<.'  multiple-step  methods  whieli  rc'cpiire  additional  intermediate  t(.'ni])oral 
data  at  the  computational  boundary.  The  implementation  ol  these  intermediate  teni])oral 
data  at  the  numerical  boundaries  lor  implicit  methods  (Til  and  T(.'l)  was  also  investigated, 
l  igure  h  depicts  a  typical  electrical  held  intensity  ol  a  tra\eling  wave  across  the  two- 
dimensional  domain  at  15  degrees  with  resj)ect  to  the  .r  coordinate  (result  ol  the  T(,'l 
Scheme).  Numerical  solutions  ol  the  three  schemes  show  conij)arable  leatures  and  are  nearly 
identical.  In  comparison  with  the  hrst-order  numerical  j>roceclure'.  no  visible  dissipatixe 
<;rrors  are  obs<;r\ecl  along  the  wave  crests. 

In  l  igure  10,  the  L2  norms  ol  calculated  electrical  lield  intensity  ol  the  SUE.  Til.  and  K.'l 
schemes  are  gi\en.  Ihe  nia.\imum  di.scTej)ancy  ol  the  L2  norm  from  the  e.xact  solution 
belongs  to  the  SUE  .scheme  with  a  magnitude  of  0.011.  The  pattern  ol  <piasi-phv.sical 
errors  exhibits  a  symmetric  structure  with  respect  to  the  i)rin.cipal  axis  ol  wave  motion 
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au<l  all  i'lrui'.s  an;  <  lu^l('^<'<l  a<lja<.<;ul  to  lh<;  uuiiK'iK  al  i)oiui<.laiy  (l  iy  lUa).  liabictilly.  tli<; 
>[)<;(- ilual ion  ol  th<'  <'xa<  l  Nolutioii  iil  iIk;  boiiii<laiy  aj>i«;ais  aN  a  p<'il inlaalioii  aiul  <-loo  not 
ut'ta'ssariiy  Ix'luii'a  to  tli<;  <;la^^  ol  initial  <l<ita  tor  wliiah  tin;  <.lill<'n;ntial  j)robl<'m  i^  \v<;ll-|)0!^<;<i. 
rh<;  .''oliition  \  ia  Hi  .s<  Ik.'iik'  yiolfl.s  a  iiiaxuiiuiii  L'J  iioiiii  lower  l>y  oin;  litill  to  that  ol  the 
SI  h  iiK'tho*]  1 1'  i”  I  he  atruetun;  ol  luiiiR'rirxil  error  i.':;  aisyniinetrie  to  the  prineipai  axiri 

ol  th«'  wa\e  [inipayation.  lloW'<'\e;.  a  <leliniti\<;  attribution  ol  this  peeuliar  tealure  to  either 
the  artilieial  nunierieal  bounUary  eoiulitions  or  the  nuiiierieal  sweep  bias  is  not  eertain  at 
pres<;nt.  lh<'  nniuer’eal  result  y,<ni<;rate<.l  by  the  I'C'l  sehenie  eontains  the  lowest  error  with 
the  iiuixinium  L‘J  norm  (U.(J02S)  just  a  lilth  ol  the  SI  E  method,  ihe  distribution  ol  quasi- 
physieal  inror  is  similar  to  tin;  two  trajiezoidal  implieit  methorls  as  exjieeted  (higs.  10b. 
lUe).  .Xomiiudly.  all  numerical  algorithms  eonsidererl  are  second-orrler  in  time  and  space. 
l.)Ut  tin;  sul.'tle  rlillen.'nce  in  num<;rical  a[)pro\imalions  lor  the  K.'l  sclieme  has  gained  it  a 
supinior  accuracy  oxer  the  ill  scheme  in  the  two-dimensional  calculation, 
lln;  eliects  ol  internn.'diate  t<'mporal  data  implementation  lor  the  apjuoximate  lactored 
sclnnin;  on  solution  accuracy  w<;n;  isolaterl  by  carrying  out  two  otherwise  idcuilical  calcu¬ 
lations  with  dillerent  data  .sets  lor  tin'  intermediate  step.  A  right-running  trans\er,se  waxe 
incliin.’rl  at  00  <legre<'s  to  tin'  .r-a.xis  was  simulated  by  the  ICl  scheme.  Numerical  results 
w<'r<'  r<'cox<;r<'<l  Isy  rx'ix.'ating  tin;  in'w  time  level  data  at  both  the  lirst  and  the  linal  nu- 
nn.'rical  sweeps  (h'ig  ila).  ainl  ijy  ch-rixiug  the  lirst  sweep  data  Irom  the  linal  time  lexel 
as  iir.licated  in  our  earlier  <,liscu.ssion  (Fig  lib),  'lhe.se  results  demonstrate  that  applying 
repetiti\<;ly  tin;  next  time  le\<’l  data  lor  l.)oth  lactored  sweeps  has  introduced  greater  errors 
than  has  applying  the  derixed  data  near  computational  bouiularies.  The  lornier  jxrocedure 
is  commonly  used  in  most  computational  Iluid  dynamics  practices  even  lor  unsteady  Iluid 
dynamics  phenomena*"*.  Its  impact  on  numerical  stability  and  lidelily  to  physics  in  .solving 
a  tinn;-dej)endent  waxe  jxrobh'in  may  be  crucial.  In  addition,  the  intermediate  solution 
Irom  an  approximate  lactored  schenn'  is  known  to  have  no  corresponding  jchysical  meaning 
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w  lii<  h  <  Im'  ail  iiii<l<'Mral>l<'  l<  atni<'  ul  ilic  jiioaii*  niij»iicit  jn<.'‘ <.'(llll<■^^.  A  loi 

ll:i^  ■■'lii.a  l<a'iiiiiiy.  il  ii<'<<»ai>.  iua\  i<'<iuii<-  I  iii‘-i  iIk-  iiiijilicU  alyonlliui  to  the 

uiiyiiuil  alii'i  iiaUU'.;  <lii'  <  l  ioii  iini'h'  ii  t<>t  mui.ii  i< 

1  hr  yiuiqi  <<1  [mai 'h'lii^  aii<tl\/<'<l  i:-  llir  <-lr<;t luinayiH'lir  waw  iinluta'cl  lay  a 

IKiiiil  <  iina'iil  SiiK  *'  <111  iMaUiapir  point  aotii<<'  will  ,^<'n<'i<itr  a  cimilai  ju<jj)a,y<itioii 

IKitW'iii  <iii'l  iIk'  i<'1<iI  i\ I'ly  .s|)<u>c  lar.^li  point  <li>tii!iu'.ion  at  tlir  iiuinl  i:otn<'<'  1.-  uimlalt' 
to  piaa\ hale  .''Ulliei<‘nl  nnnieiieal  loolntioii"* only  the  a.bility  ol  solntioni'  to  !<.‘t<un  tla.’ 
Mni|>l<'  i'a<li<t' ill.'-!  w<ui'  patta'i  n  i.-a  in\<‘.''ti“itte<l.  1  he  hoinniary  eoiialition^'  <an  tin;  ti  une<tle<.l 
<  onipMl<itioiuil  <l<aiiniin  re<[nii'e  iin  oininii  a'i.yenv'i'etoin  (o!i<’-<hiii<;nbion<il  ehaiael<;iiblic>)  laa 
\<uii>h  in  all  < lirea  t ions.  1  h<' on<'-< liiua'iisKjiuil  <  luua<;l<'ristie  l.>onn<laiy  eoiniition  is  the  e.xa'  t 
tiinl  well-paas(;<l  eoiKlilioii  lor  oiie-aliitieiisio  al  wa\e  iiiutiun.  When  apjilk'cl  to  niultiph 
(linn'iisional  [ilK'noiin'ini.  it  iiuiy  <l<‘!n'n<;i<ite  t<a  only  liist-oralcr  ac<;nt'a<;y.  (.l<;]>»;n<,ling  on  tin' 
oiienlatia'n  anal  inlrin.''ie  strna  l  nia- ol  th<;  wa\<-  iiioi iani'  ‘  ‘'h  In  j)r<.’s<'nt  an<il.\  si.s.  the 
elntraa  ta.'t ist ie  Inaiinlars'  'onalitioii.  whi>;h  's  an  inherent  an«l  inta.'iiala'al  a'ons<;<jna;iKe  aal  th<' 
llM.\-s[)lii  I  iiiy  laarnitilat  ion.  wats  iiiiposa'al  on  a'a<;h  a  ooralinata;  alirax  tion.  In  this  sa'iisa;.  this 
<ipl)roii<'h  i<'pi<’sents  th<'  worst  [laissilah'  s<<;nario  lor  .soh’.n.g  tin;  point  soiirea.'  prohia'in  in 
wliia.'h  lh<'  lai'liatinu  wa\<;  inotiaaii  is  ..i<.’seril.«'al  lay  a  (.'arta‘si«n  syst<;ui.  I  nlortunalahy.  nsin.ti, 
the  eaiaiialinaie  1 1  .uisloi iiui t ioii  to  aliyn  with  tin;  prineipal  a.\is  ol  wave  motion  has  taa  wait 
lair  I'nrtha'r  <1<'\ claapinent.  .\  t>'i)ieal  nnnieria  al  ra;snlt  using  the  T(.'l  seh<;nie  is  <!epiel<;al  in 
I'  igura’  I'J.  Ihe  :  eoiuponent  ol  the  <;l(;etrieal.  anal  ,r.  //  aonij)ona;nts  ol  tha'  luagna.’tie  ha;lal 
inla’iisil la;s  wa'i'a;  [irojeeteal  am  s<’parat<;al  plati<;s  hut  uiaintaiua;al  tit  th<;  iala.'utieai  physieal 
loe<ttion.  hiinea;  ap.i.'ntilieation  of  tha;.s<;  aoinput<;al  r<;sults  was  not  atteni])teal.  only  g<;na-ral 
leatura.’s  of  the  results  will  be  brielly  mention';  1.  At  the  singular  ])oint.  the  eurra'iil  inalu<<-s 
<i  paalariza'il  a'leetriaal  wa\(;  ni<a\ing  malialh  with  a  u<;arly  symma'lrieal  front,  wliereas  tin; 
twaa  a>amj)ouents  'af  the  magna.'tie  lia'hl  reva;al  a  laasiaally  t v  o-aliina;nsional  strueture  oj)j)osita; 
in  phasi'  across  the  point  sourca'.  i  ha’  numerical  r<;sults  r.r<  similar  to  those  ol  Analerson  s 


t  cilculalioii '  aii<i  iha  tcaluia.^  ol 

( 'alculal«.'<l  chx  liical  li<'l(l  iiitciisilio  In’  lli<'  .SI  K.  111.  anti  iC'l  .s<;lK'nlt'^  aiv  in 

l  i,t!,mt'  l-l.  .V  inuiiut hi'uniat if  puinl  suiiRc  i^  tlcsimialri.l  al  llin  inidpoijit  ol  tlie  tunipula- 
tion<il  tloniain  (i.j  —  lid)  willi  a  noiiualiztxl  wave  lon,t!,lli  ol  'Iz  anti  i«  .swittliotl  on  al  ihf 
hlarl  ol  IIk'  iiunK'iital  .sinnilation.  A  total  ol  Ui  linif  .‘>l('p.s  wore  lakon  to  allow  Uio  liis^l 
olctUital  wa\<'  lo  pa^^s  lluou,i!,h  lln-  Uuntalntl  toinpulalional  tloniain.  1  lificloif.  all  givon 
I’oull.s  an-  t'<ij)l  mt'tl  whon  the  I’chilivt.'ly  r><-n;siuvt‘  Irougli  of  tlif  wavo  Iront  inipaclt^  uj)on  tlio 
nunH'i’ital  l)oun<laiio.s.  Lila-  all  proviou^;  ono-tliint-n^^ioiuil  talfiilalion^i.  cat'll  wa\c  jiafla.’t  i.-^ 
K'xolvctl  In  nil  notk's.  iu  <  oniparc  llic  ihrcc  nuincrital  ic.-sull.s  on  an  ctjual  looting,  itlcnlical 
fonloui'  lc\<'l.''  were  u.‘^<'<l  lor  the  ii^otlynainie  plots.  In  general,  all  nuinerieal  simulations 
r<'\eal  similar  wa\’e  si  rut  lure  iu  tliat  all  isotlynamies  lorm  nearly  eoneentrie  eireles.  Small 
(lillerenees  in  the  wa\e  structure  are  probalily  intlueetl  by  relative  phase  errors  ol  each  in- 
tli\  itlual  soU’ing  stheni<.-.  luit  the  most  signilkant  tlistortion  ol  wave  sha])e  is  reeogni/etl  to 
lie  eausetl  by  lla-  oiK.’-tlimensional  charaeleristie  bountlary  eontlitions.  liie  SI.  E  inethotl. 
iK'in.g  a  single-step  e.\[>litil  im-lhotl.  protlucetl  a  mueli  better  tlelinilion  ol  the  electrical  jnil- 
salion  at  lh<'  point  source.  Hh.'  s[)ecilic  Irelun  ior  clo.sely  resenililes  the  numerical  result  ol 
,\ntl<-rson  on  a  -31)  ■  3U  nn.-sh  system'.  .Ml  implicit  schemes  tested  here  intlicaletl  a  neetl  lor 
mesh  relinemenl  at  the  singular  jioint.  However,  the  'K.'l  scheme  .seems  to  ])roduce  a  le.ss 
wave  tlistortion  than  ill  imim.-tliatels  atljacc-nt  to  the  truncatetl  com])Utational  boundaries. 
In  ortb.-r  to  tlemonslrate  the  greater  numerical  elliciency  ol  the  'K.’l  scheme,  this  proce- 
dun,'  alone  was  applh.’d  to  a  pulsating  jioint  source  inoblein  al  a  (.'LL  value  ol  2.  which  is 
unatlainal.ih'  l.n  Sl.'E  anti  ill  methotls.  The  point  source  was  moved  to  a  corner  ol  the 
calculation  tlomain  (lU.lU)  and  the  trans\er.se  wave  jracket  was  still  delined  by  2d  notles 
in  ijoth  coortlinates.  1  he  numerical  pretliction  of  the  electrical  held  intensity  is  displayed 
alP.-r  an  arbitrarily  selectetl  120  time  stejis.  The  resolution  at  the  ])oint  .source  a])])ears  to 
markedly  impro\e  but  really  is  just  highlighted  by  the  denser  contour  kwels  (Tig  11).  tV 


<[uanliili\'<.'  with  oilier  known  roulls  luiis  not  k<;<“n  achieverl  at  jnvbenl ^  Hk.' 

only  conclusion  one  may  olk'r  is  that  the  nunierical  results  are  consistent  with  the  previous 
calculations  rlepictcxl  in  k  igure  i;3.  The  accumulate*!  distortion  of  the  radiating  wave  struc¬ 
ture  generated  by  the  one-<lim*'nsionai  charact<‘ris' ic  lioundary  conditions  becomes  more 
anno.vingly  pronounced  on  r*'ctiliuear  coordinates.  ]»articularly.  in  view  ol  the  lad  that  the 
inesent  lormulation  is  designed  to  satisly  tin-  exact  one-dimensional  characteristics  ami  is 
capable  ot  aligning  the  cooixlinates  with  tlie  (uincipal  axis  ol  wave  motion,  thus  reduc¬ 
ing  the  larlield  to  a  one-<limensional  luolrlem.  ihis  obseiuation  reinforces  the  urgency  to 
complete  the  incorporation  ol  cooixlinale  translormalion  into  the  new  method. 


X 


Figure  i  l;  Ibodyiminies  of  Electric  Field  Intensity  (CFL=1.0.  TCI  Scheine.  p=4-) 


000 


5.  Conclusion 


1 


.\<'\v  trajxvuitlal  \viu<l\\ai<l  iiuiilicil  .scliciiics  oh  liic  Ilux  bplitliiig  coiicepl  were  sue- 

ee»lully  <.leveloi)e<l  to  solve  llie  lime-domain  Maxwell  e<iiialjoii.s.  Ihe  trapezoidal  consistent 
implicit  scheme  is  uncoinlitionally  stable  lor  the  linear  initial  \alue  system  and  juovirles 
a  greater  Ilexibilily  ainl  nmnerical  eflich-iiey  lo  numerically  simulate  electromagnetic  wave 
phenomena.  Applie<l  to  t  wo-<.limensional  problems.  th<'  j)res<’nt  im])licit  schemes  generated 
solutions  with  an  accurac>'  comparable  to  the  estaljlished  e.xj)licit  method  and  allowed  a 
greater  time  step  size  lo  descril.»e  the  lime-dependent  wave  phenomena. 

I  Ik.'  lull  pol<.'nlial  of  I  Ik.'  pn.'senl  eigenvector  lormulalion  is  still  not  I’ully  e.\i)lored  until  the 
general  coordinate  Iransroimation  is  I'unclional.  By  realigmueul  ol’  the  coordinate  to  the 
princi[)al  wave  motion,  the  eigenvector  rormulalion  will  be  enlumced  by  the  exact  anti  well- 
j)o.sed  oue-vlimensional.  no-rellecting  characteristic  conditions  on  truncated  coiiijvutational 
boundaries.  In  principle.  iIk-  exl<'nsion  of  the  [uocedure  to  three-dimensional  systems  is 
possible  but  may  not  b<.’  straightforward.  Ell'ort  is  ie<iuired  to  convert  the  pre.senl  apjnox- 
iiuate  factored  scheme  to  the  orginal  ADI  method  for  teni])oral  accuracy. 
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A.B.(  '  <,ucllki<.'ul  lUiiliico  ol  the  Maxwell  o<iuatioii.s 
( ’1'  L  (  'uuiaiit-r'ik'Uiiclis-Levy  uuiiilxT 

D  luaU'ix  of  the  eigen valiK' 

E  elecU'ie  liekl  intensity 

/  ,(/'  Ilux  vectors  ul  the  Maxwell  e<inations 
11  magnetic  liekl  intensity 

t.j  imlexes  ol  the  cliscretizcxl  sj)alial  points 

k  matrix  of  similar  tianslormatioii 

I  tiuK' 

I  '  <.le[)en<,lent  \arialjles  of  the  Maxwell  e<iuations 

II  (.'igen  vectors 


x,\vz  (.'artesiuu  coorcUnales 

(.  eh'ctric  permittivity 

A  eigenvalue  of  coeilicient  matrices 

//  perinea  Ijility 

E  v-C  transformed  coordinate  system 

sj  -S!/-s:  metrics  of  coordinate  transformation 


sui)erscrii)ts 

+  N'ariable  associated  with  positive  or  negative  eigenvalue  respectively 

*  Intermediate  variable  rluring  numerical  sweej) 

n  time  level 
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